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Relativistic mean field calculations have been performed to obtain nuclear density pro- 
file. Microscopic interactions have been folded with the calculated densities of finite 
nuclei to obtain a semi-microscopic potential. Life time values for the emission of proton, 
alpha particles and complex clusters have been calculated in the WKB approach assum- 
ing a tunneling process through the potential barrier. Elastic scattering cross sections 
have been estimated for proton-nucleus scattering in light neutron rich nuclei. Low en- 
ergy proton reactions have been studied and their astrophysical implications have been 
discussed. The success of the semi-microscopic potentials obtained in the folding model 
with RMF densities in explaining nuclear decays and reactions has been emphasized. 



I. INTRODUCTION 

Relativistic mean field (RMF) approach has 
proved to be very successful in explaining 
different features of stable and exotic nuclei 
such as ground state binding energy, deforma- 
tion, radius, excited states, spin-orbit split- 
ting, neutron halo, etc[l[. Particularly, the 
radius and the nuclear density are known to 
be well reproduced in this model. In nuclei far 
away from the stability valley, the single par- 
ticle level structure undergoes certain changes 
in which the spin-orbit splitting plays an im- 
portant role. Being based on the Dirac La- 
grangian density, and thus naturally incorpo- 
rating the spin degree of freedom, RMF is par- 
ticularly suited to investigate these nuclei. 

RMF calculations have been found to pro- 
vide good description of densities, particu- 
larly at large radii. Processes, such as pro- 
ton and alpha decay, low energy scattering 
and reaction, probe this part of the nuclear 
volume. Thus, a good description of such 
processes may be expected if we use micro- 
scopic nucleon-nucleon (NN) interactions and 
densities from RMF calculations to construct 
nuclcon(nucleus)-nucleus potentials. This will 
be very useful to extend any calculation to ar- 
eas far from the stability valley, where data 
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are scarce and difficult to obtain in near fu- 
ture. For example, low energy reactions are 
very important from the astrophysical point of 
view. In astrophysical environments, neutron 
and proton reactions are the keys to nucle- 
osynthesis of heavy elements. However, very 
often, the target nuclei are not available in ter- 
restrial laboratories and one needs to depend 
on theoretical inputs. 

The presentation has been arranged in the 
following manner. In Section 2, the method 
of constructing the potential is outlined. The 
next section presents a few theoretical results 
for density, followed by selected results of cal- 
culation for decays involving emission of pro- 
ton, alpha and complex clusters. We also 
briefly discuss the results of elastic proton 
scattering in neutron rich nuclei. A few low 
energy proton reactions in A = 60 — 80 region 
have also been studied and astrophysical im- 
portance of such reactions has been discussed, 
Finally, we summarize our results. 



II. METHOD 

There is a large body of work for the top- 
ics that have been presented in this talk. In 
this section, we restrict ourselves to the salient 
points of the method of calculation that have 
been followed to obtain the semi-microscopic 
potentials used in the present work. 

Theoretical density profiles have been ex- 



tracted from RMF calculations. In this ap- 
proach, nuclcons interact via exchange of a 
number of mesons. There arc different vari- 
ations of the Lagrangian density as well as 
a number of different parametrizations. In 
the present work we have employed the FSU 
Gold[2| Lagrangian density. It contains, apart 
from the usual component describing a sys- 
tem of nucleons interacting via exchange of 
mesons, nonlinear terms involving self cou- 
pling of the scalar-isoscalar and the vector- 
isoscalar meson, as well as coupling between 
the vector-isoscalar meson and the vector- 
isovector meson. 

In the conventional RMF+BCS approach 
for even-even nuclei, the Euler-Lagrange equa- 
tions are solved under the assumptions of clas- 
sical meson fields, time reversal symmetry, no- 
sea contribution, etc. Pairing is introduced 
under the BCS approximation. Since accuracy 
of the nuclear density is very important in our 
calculation, we have solved the equations in 
co-ordinate space. The strength of the zero 
range pairing force is taken as 300 MeV-fm 
for both protons and neutrons. These values 
have been chosen to represent a good fit for 
the binding energy values. In nuclei contain- 
ing odd number of neutrons or protons, the 
tagging approximation has been used to spec- 
ify the level occupied by the last odd nucleon 
of either type. We have observed that mod- 
erate variations of the pairing strength do not 
influence the results to any great extent. 

Effective NN interactions, such as density 
dependent M3Y (DDM3Y) |343 , or that of 
Jcukenne, Lejeune and Mahaux (JLM)@ may 
be used to construct the nucleon (nucleus)- 
nucleus potential. Both the interactions men- 
tioned above have been derived from nuclear 
matter calculation and have been applied in 
finite nuclei with success. 

The DDM3Y interaction [1 Q is obtained 
from a finite range energy independent M3Y 
interaction by adding a zero range energy de- 
pendent pscudopotcntial and introducing a 
density dependent factor. The density depen- 
dence may be chosen as exponential[3j or be of 
the form C(l - f3p 2 / 3 )\^. The constants were 
obtained from nuclear matter calculation [fl 



as C = 2.07 and /3 = 1.624 fm 2 . We have 
adopted all the standard parameters in our 
calculation without any modification. This 
interaction has been employed widely in the 
study of nucleon nucleus as well as nucleus 
nucleus scattering and radioactivity. 

In the JLM potential 0, finite range of the 
interaction has been introduced by including 
a Gaussian form factor^- Here, we have 
adopted the global parameters for the inter- 
action and the default normalizations 0] ■ 

The semi-microscopic nuclear potentials 
have been obtained by folding the interac- 
tions with the microscopic densities obtained 
in the RMF calculation. The Coulomb po- 
tential has similarly been obtained by folding 
the Coulomb interaction with the microscopic 
proton densities. The spin-orbit potential was 
chosen from the Scheerbaum prescription Q. 
The total potential consists of the nuclear 
part, the Coulomb potential as well as the 
centrifugal potential. We have not included 
the contribution of isovector (Lane) potential. 
However, we expect its effect to be small. 

III. RESULTS 

A. Density 

The charge density has been obtained from 
the point proton density p p by convoluting 
with a Gaussian form factor to account the 
finite size of the proton. We plot in Fig. [T]the 
calculated charge density for 62 Ni and 66 Zn 
as representatives of our results. One can see 
that the theoretical and experimental values 
agree very well, particularly at larger radii val- 
ues, which is the region expected to contribute 
to the optical potential at low energy. Other 
nuclei also show similar agreement. 

B. Decay 

Life time values have been calculated for dif- 
ferent types of decays, i.e. proton, alpha and 
cluster radioactivity. We employ the super- 
asymmetric fission model (SAFM). The decay 
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FIG. 1: Calculated charge density in Ni and 
66 Zn (solid lines) compared with experimental 
measurements (filled circles). 



product is assumed to tunnel through the bar- 
rier created by the repulsive Coulomb (and the 
centrifugal) potential and attractive nuclear 
potential, constructed by folding. The prob- 
ability of barrier penetration has been calcu- 
lated in the WKB approximation. The assault 
frequency is obtained from the zero-point vi- 
bration energy Q that, in turn, has been cal- 
culated from the Q-values. Unlike some other 
works, we have not normalized the optical po- 
tential, and have used the values that were 
obtained from nuclear matter calculations, for 
the decay studies. 

Many calculations exist for proton 
radioactivity [uj- Our results for the proton 
radioactivity half life calculation, in nuclei 
where such decay has been observed from 
low energy states, are tabulated in Tabic Q] 
For comparison, we have also presented the 
results for the DDM3Y interaction of Basu et 
al.jllj. who have used a simple phenomeno- 
logical form of the densities. We have also 
shown the uncertainties in the calculated half 
life values corresponding to the errors in the 
measured Q-values within parentheses. Our 
results deviate from experimental measure- 
ments in 147 Tm, 150 Lu and 156 Ta, and more 
prominently in 185 Bi and 177 T1*, the last two 
results being off by an order of magnitude. 
One can show that the last two discrepancies 
can be explained as the effect of configuration 
mixing [1 Of . 

It was often suggested that the half life 
values for proton radioactivity of 109 I and 



TABLE I: Proton decay half lives (T) for spher- 
ical proton emitters with even neutron number. 
Results are for the DDM3Y interaction. The an- 
gular momentum of the proton involved is given 
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H2,ii3q s cou 1q; no i he reproduced without the 
inclusion of deformation effects. For example, 
a deformation of the order of f3 ~ 0.05 — 0.15 
was judged to be essential to explain the ob- 
served data. However, our calculation repro- 
duces the data for 109 I with considerable accu- 
racy. In the deformed calculations, it is usu- 
ally assumed that the deformation of the par- 
ent and the daughter nuclei are identical. The 
daughter in this particular case is 108 Tc. Tc 
nuclei are well known vibrational nuclei with 
very small deformation. One possibility may 
be that the deficiency of the Woods Saxon po- 
tential far away from the stability valley was 



responsible for the failure of the earlier cal- 
culations. We also stress that our results are 
nearly identical for both the NN interactions. 
The results for 113 Cs and to some extent 112 Cs 
are not reproduced so well, which may be an 
effect of deformation. However, in none of 
them do we have an order of magnitude dis- 
agreement between theory and experiment as 
obtained in certain other calculations. 

Alpha decay have been studied in Refs. fl2l - 
[lil ]. For emission of complex products, there is 
another factor, the spectroscopic factor, intro- 
duced to incorporate the preformation prob- 
ability. It contains the nuclear structure ef- 
fects, and may be thought as the overlap be- 
tween the actual ground state configuration 
of the parent and the configuration described 
by the complex decay product coupled to the 
ground state of the daughter. Obviously, it is 
expected to be much less than unity as there 
are contributions from many other configura- 
tions other than the one mentioned above. 

In a small mass region, we do not expect 
the spectroscopic factor to vary to any large 
extent. For example, in the superheavy nuclei 
presented in Fig. [2l we have taken a constant 
value 1.4 x 10 -2 for all the decays from a fit of 
the half life values [HI]. In some other calcula- 
tions, the spectroscopic factors have been cal- 
culated from theory or phenomenological for- 
mulas have been obtained for them[l3f. 



study emission of complex clusters [H, EH- 
Basu 171 has studied the nuclear cluster ra- 



dioactivity also in the framework of SAFM us- 
ing a phenomenological density and the real- 
istic M3Y interaction. Bhagwat and Gamb- 
hir [ll| have used densities from RMF calcu- 
lations. It has been suggested fl§j that in the 
case of decay of heavy clusters, the spectro- 
scopic factor may scale as 



S = (S Q ) 



(A-l)/3 



(1) 



where A is the mass of the heavy cluster and 
S a is the spectroscopic factor for the a-decay. 
Thus a plot of log 10 S against A should be 
a straight line. In Fig. [31 we have plotted 
the negative of log 10 S for the decays where 
both the parent and the daughter are even- 
even nuclei against the mass number of the 
cluster and plotted a best fit line. One can 
see that the points fall nearly on a straight 
line. 
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FIG. 2: Calculated and experimental half life val- 
ues in superheavy nuclei 

The above method may be extended to 
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FIG. 3: Negative of logarithm of spectroscopic 
factors (S) as a function of cluster mass number 
A for even-even parents and daughters. 



C. Elastic scattering 

Elastic proton scattering in inverse kine- 
matics provides a test for the calculated den- 
sities, particularly in absence of electron scat- 
tering data. Density information is not avail- 
able in exotic neutron rich nuclei in very light 



mass region. Elastic scattering cross sections 
for scattering of these nuclei from proton tar- 
get have been calculated[lj| with the opti- 
cal model potential generated in the semi- 
microscopic approach. In Fig. [H the mea- 
sured partial cross sections for proton scatter- 
ing of 6,8 He in inverse kinematics have been 
compared with theory. The real and the imag- 
inary parts of the potential have been chosen 
as 0.8 times and 0.2 times the folded value for 
the DDM3Y potential. One can see that the 
experimental values are nicely reproduced. 



D. Low energy proton reactions of 
astrophysical significance 

Capture and charge exchange reactions at 
very low energy play a very important role 
in nucleosynthesis. Particularly, rapid pro- 
ton capture (rp) process in explosive nucle- 
osynthesis is a basic ingredient in driving the 
abundance along the N = Z line[2(|. As this 
process has to overcome a large Coulomb bar- 
rier, it can occur only at a higher tempera- 
ture range. For example, X-ray bursts pro- 
vide large fluxes of protons at peak tempera- 
tures around 1-2 GK and are expected to play 
a significant role in the creation of nuclei up 
to mass 110. 

The rp-process proceeds along the N = Z 
line in mass 60-80 region. In nature, the 
important proton capture reactions may in- 
volve certain nuclei as targets which are not 
available to us. Hence, experimental infor- 
mation is difficult, if not impossible, to ob- 
tain, at least in near future. In such a sit- 
uation, one has to rely on theory for the re- 
action rates. Cross sections were calculated 
for proton capture reactions in mass 60-80 
region using semi-microscopic optical poten- 
tials in the local density approximation with 
phenomcnological density prescriptions. How- 
ever, far from the stability valley, these pre- 
scriptions may not represent the actual densi- 
ties very well, leading to considerable uncer- 
tainties in the reaction rates. Very often, the 
reactions rates are varied by a large factor to 
study their effects. For example, Schatz(2~i| 



varied the rates of certain reactions by a fac- 
tor of one hundred. Obviously, this makes the 
results uncertain to some extent and affects 
the final abundance. 

A fully microscopic calculation may be used 
to estimate the rates to reduce the above un- 
certainty. A consistent framework for calcu- 
lation may be constructed based on micro- 
scopic densities and may be extended to un- 
known mass regions with confidence. In the 
present work, we have tried to calculate the re- 
action rates taking densities from a purely mi- 
croscopic model, i.e. RMF. We have already 
mentioned that RMF is particularly suitable 
to describe nuclei far away from the stability 
valley where experimental knowledge is scarce. 
A semi-microscopic optical potential obtained 
by folding an appropriate microscopic NN in- 
teraction is expected to be more accurate and 
may do away with the necessity of any arbi- 
trary variation in the reaction rates. 

However, even in a semi-microscopic optical 
potential, there often remain certain parame- 
ters that can be fixed only after comparison 
with experiment. We have compared the re- 
sults for a number of (p, 7) and (p, n) reactions 
in mass region A = 60 — 80 for which exper- 
imental cross sections are available. This has 
helped us in determining a set of parameters 
for this mass region. 

The reaction calculations have been per- 
formed with the computer code TALYS 1.2[22j 
assuming spherical symmetry for the target 
nuclei. The real part of the potential has been 
obtained by normalizing the folded DDM3Y 
potential by a factor of 0.7, while the imag- 
inary part, by a factor of 0.1, so as to ex- 
plain the S-factors obtained in the above ex- 
periments. We have employed the full Hauser- 
Feshbach calculation with transmission coeffi- 
cients averaged over total angular momentum 
values and with corrections due to width fluc- 
tuations. The gamma ray strength has been 
calculated in the HFB and HF+BCS model. 

In Figs. [5] and [HI we plot the results of our 
calculations for S-factors and compare them 
with experimental values for (p, 7) reactions 
for various targets in mass 60-70 region. The 
energy relevant to the rp-process in this mass 
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FIG. 4: Partial cross section for elastic proton scattering in inverse kinematics. The projectile energies 
of 6 He and 8 He are 71 MeV/A and 72 MeV/A, respectively. 



region lies between 1.1 to 3.5 MeV. As the 
cross-section varies very rapidly at such low 
energy a comparison between theory and ex- 
periment is rather difficult. The usual practice 
in low energy nuclear reaction is to compare 
another key observable, viz. the S-factor. It is 
given by 

S(E) = Ea{E)e 27Tri (2) 

where E is the energy in centre of mass frame 
in keV, cr(E) indicates reaction cross-section 
in barn and rj is the Sommcrfcld parameter 
with 

27^ = 31.29^^/1; (3) 
V E 

Here, Z p and Z t are the charge numbers of 
the projectile and the target, respectively and 
\i is the reduced mass (in amu). The quantity 
S-factor varies much more slowly than reac- 
tion cross-sections as the exponential energy 
dependence of cross-section is not present in 
it. For this reason, we calculate this quantity 
and compare it with experimentally extracted 
values. 

In Fig. [Jj we plot the results for cross sec- 
tions to populate the ground state and the 
first two excited states for 63,65 Cu(p, 7) reac- 
tions. Results for proton capture by other nu- 
clei as well as charge exchange reactions have 
also been calculated in this mass region and 
found to agree reasonably well with experi- 
mental observations [23-123. 



Once the parameters have been fixed, we 
employ them to calculate the rates of a num- 
ber of astrophysically important proton re- 
actions. Some very important N — Z nu- 
clei, which have the highest abundance in an 
equilibrium in a chain, are termed as wait- 
ing points [2(3] for the chain. These nuclei 
have negative or small positive Q-values for 
proton capture. An equilibrium between the 
(p,7)/(7,p) processes is established and the 
rp process may have to wait for beta-decay or 
a-capture to proceed to heavier nuclei. Cer- 
tain N = Z waiting point nuclei with A < 80, 
viz. 64 Zn, 68 Se, 72 Kr, and 76 Sr have long half 
lives, their total lifetime being large compared 
to the time scale of typical X-ray bursts (10- 
100 sec). Thus, they may produce a bottle- 
neck in the rp-process that would slow down 
the rate of hydrogen burning and necessitate 
extended burst tails unless two proton cap- 
ture can reduce these half lives and bridge the 
waiting points. X-ray burst model calcula- 
tions are therefore particularly sensitive to the 
rates of proton capture for these nuclei. We 
have used the microscopic approach, outlined 
in the present work, to calculate the rates with 
an aim to study the bridging of the waiting 
point nuclei. 

For calculation of proton capture at wait- 
ing points, a small network, which includes 
the following processes, has been employed. 
The waiting point nucleus with Z = N, which 
acts as a seed, may capture a proton. The re- 




FIG. 5: S-factors extracted from experimental measurements compared with theory for ' Ni and 
63,65q u targets. Solid and dashed lines indicate respectively the results of the HF+BCS and HFB 
approaches for level density and El gamma strength. 




FIG. 6: S-factors extracted from experimental measurements compared with theory for ' Zn. See 
caption of Fig. [S] for details. 



suiting nucleus, with Z = N + 1, may either 
capture another proton or undergo photodis- 
integration emitting a proton to go back to 
the seed nucleus. The nucleus with Z = N + 2 
may also undergo photodisintegration. In ad- 
dition, all the three nuclei mentioned above 
may undergo /3-decay. The photodisintegra- 
tion rates at different temperatures have been 
calculated from the proton capture rates using 
the principle of detailed balance. The density 
has been taken as 10 6 gm/cm 3 unless other- 
wise mentioned. The proton fraction has been 
assumed to be 0.7. 

Fig. [8] shows the change in the effective half 
life of Ge in explosive hydrogen rich envi- 
ronment. For comparison, we have also plot- 
ted the results calculated from the rates in 



Rauschcr ct al. 28] by dash-dotted lines. The 
effects of the uncertainties on the half life val- 
ues in the Q-values have been indicated in the 
figures by dotted lines. The half life decreases 
and possibly goes to a value substantially less 
than ten seconds, the minimum duration of an 
X-ray burst. However, one sees that the un- 
certainty in mass measurement prevents one 
from reaching any firm conclusion. Depend- 
ing on the actual value of the masses, it may 
even be possible that a burst of the order of 
ten seconds cannot bridge this waiting point 
effectively. 

A larger network, starting from 56 Ni seed 
and consisting of (p, 7), (j,p) reactions and 
/3-decay, has been employed to study the the 
time evolution of abundance. In Fig. [S]we plot 
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FIG. 7: Partial cross sections for (p, 7) reactions 
to the low- lying states in Zn for 63 Cu (left panel) 
and 65 Cu (right panel) targets. Squares, circles 
and triangles represent data for transition to the 
ground state, the first excited state (multiplied 
by 10) and the second excited state (multiplied 
by 100), respectively. 
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FIG. 9: The rp-process path for 1.2 GK. 
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FIG. 8: Effective half life values of 64 Ge as a func- 
tion of temperature. The solid line represents the 
results of our calculation while the dotted lines 
mark the two extremes for the errors in the Up- 
values of the reactions involved. The dash dot- 
ted line shows the results obtained using the rates 

from m. 



the path followed by the rp-process between 
A = 56 and A = 80 for T = 1.2 GK. The dark 
boxes indicate the waiting point nuclei. The 
elements are indicated at the left of the dia- 
gram while the neutron numbers are shown at 
the bottom. The lines indicate the path fol- 
lowed by nucleosynthesis. The paths through 
which more than 10% of the total flux flows are 



indicated by black lines while gray lines show 
the corresponding paths for flux between 1% - 
10%. One can see that 64 Ge is easily bridged 
by two proton capture at 1.2 GK. However, 
we have found that at a higher temperature, 
viz. 1.5 GK, this waiting point delays the pro- 
cess so that the /?-decay of 64 Ge contributes 
significantly. At the next waiting point, 68 Sc, 
the photodisintegration is sufficiently strong 
so that, independent of the temperature, /?- 
decay is practically the only available path. 
This delays the nucleosynthesis significantly. 
At 72 Kr, the inverse process predominates in 
higher temperature driving the flux through 
/3-decay. Thus, here also, lower temperature 
helps nucleosynthesis speed up. The waiting 
point at 76 Sr presents a different picture where 
the path essentially does not depend on the 
temperature and principally flows along de- 
cay. It is clear that the actual process is sig- 
nificantly dependent on the model of the burst 
process where the temperature and the proton 
fraction are functions of time. 

The abundance at two different tempera- 
tures are plotted as a function of time in Fig. 
[TO] We find that at 1.2 GK, at the end of 100 
seconds, the population that reaches A = 80 
or beyond is more than 1.5 times than the 
corresponding quantity at the higher temper- 
ature of 1.5 GK. In both the cases, the popu- 



lation beyond A = 76 is significant. It should 
be noted that these numbers are strongly de- 
pendent on the masses of the nuclei near the 
waiting point. Masses of these nuclei have ei- 
ther not been measured, or measured very im- 
precisely. Thus one has to depend on theory. 
In fact, the final abundance at various masses 
change by a large factor as one uses estimates 
from different formulas. This conclusion is in 
agreement with the variation in effective half 
life values seen in Fig. [8] In the present in- 
stance, the Duflo-Zuker formula[29| has been 
used. 
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FIG. 10: Evolution of abundance of mass at the 
waiting points in explosive proton-rich astrophys- 
ical environments. 



IV. SUMMARY 

The techniques used to obtain semi- 
microscopic optical potentials from theoretical 



densities using RMF and effective NN interac- 
tions have been described. Life time for emis- 
sion of protons, neutrons and complex nuclei 
from nuclei have been successfully calculated 
in the WKB approximation using the real part 
of the above potential. The potential has also 
been used to describe elastic scattering and 
proton reactions at low energy. Astrophysical 
implications of the results of proton capture 
reactions have been discussed. Representative 
results for above calculations have been pre- 
sented. It is possible to conclude that semi- 
microscopic optical potentials obtained in the 
folding model with RMF densities are very 
much successful in explaining nuclear decays 
and reactions over the entire mass region. 
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